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Abstract 

High-performance electrothermal thrusters 
operate In a low nozzle-throat Reynolds number 
regime. Under these conditions, the flow boundary 
layer occupies a large volume Inside the nozzle, 
contributing to large viscous losses. Four nozzles 
(a conical, bell, trumpet, and modified trumpet) 
and a sharp-edged orifice were evaluated over a 
Reynolds number range of 500 to 9000 with unheated 
nitrogen and hydrogen. The nozzles showed signifi- 
cant decreases In specific Impulse efficiency with 
decreasing Reynolds number. At Reynolds numbers 
less than 1000, all four nozzles were probably 
filled with a large boundary layer. The discharge 
coefficient decreased with Reynolds number In the 
same manner as the specific Impulse efficiency. 

The bell and modified trumpet nozzles had discharge 
coefficients 4 to 8 percent higher than those of 
the cone or trumpet nozzles. The Two-Dimensional 
Kinetics (TDK) nozzle analysis computer program 
was used to predict nozzle performance. The 
results were then compared to the experimental 
results In order to determine the accuracy of the 
program within this flow regime. 

Introduction 

Low Reynolds number nozzles have large flow 
boundary layers and associated viscous losses that 
decrease the performance of low- thrust electro- 
thermal propulsion devices. Thrust losses approach 
30 percent In the Reynolds number range that is 
characteristic of high-performance reslstojets and 
arcjets.! The one-dimensional Isentroplc expan- 
sion analysis used for evaluating nozzle perfor- 
mance losses generally agrees with experimental 
data within a few percent for high Reynolds number 
nozzles ( >1 0 000). Because of the effects of large 
boundary layers at low Reynolds numbers, the Isen- 
troplc expansion analysis Is no longer accurate. 

Earlier work In the area focused primarily on 
conical nozzles of varying half-angles and area 
ratios. It was shown that the boundary layer and 
associated viscous losses cause a decrease In 
thrust coefficient with decreasing Reynolds number. 
Three loss mechanisms were Identified, a boundary- 
layer effect, a heat-transfer effect, and flow- 
divergence losses. The viscous losses are usually 
much larger than divergence losses. Low Reynolds 
number nozzle flow has an Invlscld core displaced 
from the wall by a viscous boundary layer. Flow 
velocity In both regions Is expected to be less 
than Ideal. 2 Viscous and divergence losses at 
throat Reynolds numbers less than 1000 can reduce 
the overall thrust power efficiency to less than 
75 percent.! when heated nozzle flow Is eval- 


uated, radiation and conduction losses can account 
for 10 to 20 percent of the Input gas power. 3 
In vacuum chambers, where the ambient pressure Is 
sufficiently high, heat-transfer losses due to 
convection can also occur. 4 * 5 Under the condi- 
tions of large boundary layers and viscous losses, 
theoretical analysis of nozzle performance using 
the assumption of one-dimensional Isentroplc 
expansion Is no longer accurate and must be 
modified." 1 .2.6 

NASA has Initiated In-house and contracted 
efforts to develop a high-performance, storable 
propellant reslstojet. The potential for perfor- 
mance gain by minimizing nozzle losses has been 
identified as an area for further Investigation. 

The purpose of the work presented in this paper 
was to extend the data base to Include nozzles 
with varying outlet geometries and to examine only 
boundary- layer effects by using unheated gas. 

Four nozzles of varying divergent section 
geometries - conical, bell, trumpet, and modified 
trumpet - and a sharp-edged orifice were tested 
with nitrogen and hydrogen over the Reynolds number 
range of 500 to 9000. 

Present analytical methods used for nozzle 
analysis have not been validated at low Reynolds 
numbers. The Two-Dimensional Kinetics (TDK) nozzle 
analysis program was run, and Its results were 
compared to the experimental data. The TDK pro- 
gram evaluates two-dimensional and viscous effects 
on the performance of liquid propellant exhaust 
nozzles. 7 It should be noted that TDK was written 
for high Reynolds number flow and has not been 
validated In the low Reynolds number regime. The 
version of TDK used for this paper was version 2.4, 
December 1984. 

The TDK program calculates gas composition In 
the nozzle chamber and throat, assuming chemical 
equilibrium. By using this Information, the pres- 
sure profile within the chamber and throat are 
calculated from Isentroplc flow relations. Then, 
the one-dimensional nonequilibrium flow relations 
are Integrated for the flow beginning at the con- 
verging section of the nozzle and ending at an 
axial station located beyond the throat plane. 
During this calculation, the program Is Instructed 
to assume frozen chemical composition. The results 
obtained In the throat section are then used to 
estimate two-dimensional effects In the transonic 
region of the nozzle throat. The purpose of these 
calculations Is to approximate an Initial data 
line across the nozzle throat In order to start 
the method-of-characterlstlcs calculations. The 
method of characteristics Is used to determine the 
properties of the flow within the divergent nozzle 
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and thereby calculate the loss In nozzle perfor- 
mance caused by flow divergence. The flow condi- 
tions along the nozzle wall, which are calculated 
In the method of characteristics, are used to 
perform a boundary-layer analysis. From this anal- 
ysis, the decrement in performance due to viscous 
effects Is obtained and subtracted from the Invls- 
c 1 d performance to obtain the final results. 

The benefit of using TDK over simple one- 
dimensional Isentroplc flow relations Is that It 
considers gas species, divergence losses, and 
viscous losses. Viscous losses are determined by 
calculating a boundary-layer displacement thick- 
ness, assuming viscous effects at the wall, and 
subtracting this performance loss from the Invlscld 
core performance calculation. 

This paper presents specific Impulse efficien- 
cies and discharge coefficients of the four nozzles 
evaluated over a wide range of throat Reynolds 
numbers. Experimental measurements and analytical 
results are compared. 

Apparatus 

Nozzles 

Four converging-diverging nozzles with differ- 
ent diverging contours were tested along with a 
sharp-edged orifice. Nozzle 1 was a 20° half-angle 
cone with an area ratio of 120:1 and a throat diam- 
eter of 0.653 mm. Nozzle 2, a bell-shaped nozzle 
with a throat diameter of 0.711 mm, and nozzle 3, 
a trumpet with a throat diameter of 0.671 mm, had 
area ratios of 150:1 and 125:1, respectively. 

Nozzle 4 was a modified trumpet shape - the diver- 
ging section near the throat was a 15° half-angle 
cone that expanded Into a trumpet with a 135:1 
area ratio and a throat diameter of 0.640 mm 
(Fig. 1 and Table 1). The orifice was a flat 
plate with a sharp-edged 0.724-mm-dlameter hole. 

Nozzle Interchangeability was simplified by 
using a quick-connect flange assembly with a rubber 
o-rlng (Fig. 2(a)). The nozzles were supplied by 
Rocket Research Corporation In conjunction with an 
ongoing high-performance reslstojet contract. 

Thrust Measurement 

Thrust measurements were obtained using a 
thrust stand capable of measuring thrust levels as 
low as 2.9 mN. The thrust stand consisted of a 
horizontal mounting plate supported by a double 
set of four flexure plates that allowed motion In 
the horizontal direction (Fig. 2(a) and (b)). The 
flexures were designed so that a small force pro- 
duced a large displacement. This helped reduce 
two major thrust-measurement problems - thermal 
drift and vibration. The thermal expansion of the 
thrust stand components was small compared to the 
movement of the thrust stand for low-thrust tests, 
and the flexures had a low resonant frequency, 
preventing any high-frequency vibrations from 
reaching the thruster mounting plate. The low- 
frequency vibrations were removed by use of a mag- 
netic damper assembly. The damper consisted of a 
permanent magnet physically attached to the 
thruster mounting plate and inserted Into an 
annular electromagnet. The output of the electro- 
magnet was dependent on the movement of the 
thruster mounting plate. Thrust plate movement 
was detected by monitoring the output of the 
linear variable displacement transducer ( LVDT) . 


When movement was detected, the LVDT output was 
amplified and used to energize the electromagnet. 

A magnetic field was produced which applied a force 
on the permanent magnet opposing the thrust plate 
motion. 

Various thrust ranges were measured by varying 
the stiffness of the flexures and, If required, 
the stiffness of the propellant feed tube. The 
feed tube, a 3. 2-mm-dlameter thin-walled stainless 
steel tube, was brought onto the thrust stand per- 
pendicular to Its axis of movement, thus acting as 
an additional flexure. 

Calibration was performed by adding known 
weights to a pulley assembly that deflected the 
thrust stand a given amount. Thrust stand tares 
were small and reproducible, and the estimated 
thrust measurement error was less than 5 percent 
at low thrust (<66 mN) and less than 2 percent at 
higher thrust levels. 

Windage effects, the Influence of circulating 
gases In the test chamber on thrust measurement, 
were reduced to less than 0.1 mN by Installation 
of a windshield between the nozzle and the thrust 
stand body. These effects were measured by flowing 
gas through an orifice that was mounted close to, 
but not attached to, the horizontal thruster mount- 
ing plate, and checking for any thrust plate 
deflection. 

Vacuum Chamber 

The vacuum chamber used had a diameter of 
1.0 m and a length of 2.0 m. A 30-cm pipe con- 
nected the chamber to a 1x10 s llter/mln lobe rotary 
blower that was backed by a 2xl0 4 llter/mln rough 
pump. Typical vacuum chamber pressures varied from 
1.06x10-1 to 66.5 N/m 2 (8xl0‘ 4 to 0.5 torr) 
depending on the thrust level selected. Tests with 
low Reynolds numbers ( <2000) had vacuum chamber 
pressures of less than 1.33 N/m 2 (10~ 2 torr). 
Pressure ratios, nozzle plenum pressure divided by 
vacuum chamber pressure, were kept above 10 3 to 
prevent the possibility of shock formation In the 
nozzle. 8 All measured thrust values were cor- 
rected for background pressure which amounted to 
less than 1 percent of the measured thrust. 

Instrumentation 

Gas flow was measured by using laminar-flow, 
heated-tube flowmeters that are thermally sensitive 
to flow rate. For a constant tube reference tem- 
perature, the gas flow rate is proportional to the 
current Input to the tube. The low range, 0 to 5 
standard llter/mln, flowmeter was calibrated using 
a soap-bubble calibrator. 9 The high-range, 0 to 
50 standard llter/mln, flowmeter was calibrated by 
filling a known volume to a certain pressure over 
a measured time Increment. Errors for the low- 
flow-rate meter were less than 1 percent, and, for 
the high-flow-rate meter, less than 3 percent. 

Pressures were measured by a stainless steel 
diaphragm transducer. The pressure tap was located 
approximately 12-mm upstream of the nozzle throat. 
The calibration was checked before every test, and 
errors In measurement were less than 1 percent. 
Vacuum-chamber pressures were measured with a cold 
cathode-ion gauge; estimated measurement errors 
were +20 percent for low pressure (<0.133 N/m 2 
(10‘ 3 torr)). Temperatures were measured by a 
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stagnation probe with a chromel -al umel thermo- 
couple mounted In the gas stream approximately 
12-mm upstream of the nozzle throat; estimated 
measurement errors were less than ±3 °C. 

Procedure 

The test matrix consisted of three variables - 
throat Reynolds number, nozzle contour, and propel- 
lant. The Reynolds number Is based on throat 
conditions and Is equal to 4 times the mass-flow 
rate divided by the product of it, viscosity, and 
throat diameter. There were four nozzles - a coni- 
cal, bell, trumpet, and modified trumpet - and a 
sharp-edged orifice. The propellants used were 
unheated nitrogen and hydrogen. Since the tests 
were run at constant conditions, the Reynolds 
number was varied by changing the mass-flow rate. 
The Reynolds number range tested, 500 to 9000, 
spanned two thrust ranges, 2.8 to 12.0 mN and 12.0 
to 60 mN. 

The thrust stand was set up for one thrust 
range, and a propellant was selected. The nozzle 
was attached to the thrust stand assembly. Its 
exit was plugged, and the fittings were checked 
for leaks. After ensuring that the system was 
leak tight, the vacuum chamber was evacuated. A 
mass-flow rate was set with a thrust zero taken 
before and after each thrust measurement. Thrust 
measurements were recorded twice per Reynolds 
number set point (flow rate). Each nozzle was run 
three times In random order to accurately define 
the thrust measurement deviations, and a thrust 
stand calibration check was performed before and 
after each nozzle change. 

This procedure was followed for both propel- 
lants, over both thrust ranges, and for all noz- 
zles. The quantities recorded were thrust, Inlet 
gas pressure, Inlet gas temperature, test cell 
pressure, and mass-flow rate. Thrust was cor- 
rected for test cell pressure using 

h = F m + P a A e (1) 

where F c Is the corrected thrust, F m Is the 
measured thrust, P a Is the ambient test cell 
pressure, and A e Is the nozzle exit area. 

Thrust corrections were less than 1 percent over 
the entire range. 

The specific Impulse and thrust coefficient 
were determined using F c . Specific impulse 
efficiency Njsp 1$ the ratio of actual specific 
impulse (calculated from measured thrust and mass- 
flow measurements) to theoretical maximum specific 
Impulse. Theoretical maximum values of specific 
Impulse at room temperature for hydrogen and nitro- 
gen are 300 and 78 sec, respectively. The specific 
Impulse efficiency measures the extent of the 
velocity losses due to nozzle viscous and divergent 
effects. The discharge coefficient Cq Is the 
ratio of measured mass-flow rate to theoretical 
maximum mass-flow rate. For sonic conditions at 
the nozzle throat, Cq Is dependent on the 
throat boundary-layer thickness. The discharge 
coefficient is therefore a measure of the loss of 
mass-flow due to viscous effects from the nozzle 
throat region. For a hydrogen or nitrogen propel- 
lant at room temperature, the thrust coefficient 
can be expressed as 

C F = 1.7 C D N Is p (2) 


for nozzle area ratios between 100 and 200. 6 
Furthermore, for an orifice plate, Eq. (2) also 
applies. Variations In the thrust coefficient can 
be readily obtained from figures displaying Nx$p 
and Cq as a function of Reynolds number (Figs. 3 
to 6). 

The assumptions used for the TDK analysis are 

(1) There Is frozen chemical composition. 

(2) There Is no loss of mass from the system. 

(3) Each component of the gas Is a perfect gas. 

(4) The flow Is axlsymmetrlc (In the nozzle throat 
region the flow Is dependent on the local wall 
geometry only) . 

(5) There are wall viscous effects. 

(6) Core flow calculations are performed assuming 
Invlscld, compressible gas. 7 

Discussion of Results 

The major loss mechanism In low Reynolds num- 
ber nozzles Is viscous flow loss associated with 
the flow boundary layer, which occupies a large 
volume Inside the nozzle. Data exploring these 
losses are limited, particularly In cases of noz- 
zles with different divergent contours. Most of 
the past work was performed with heated gas, which 
can Introduce an additional heat-transfer loss 
mechanism. The heated-gas tests represent a more 
accurate comparison with actual low-thrust, high- 
performance reslstojets and arcjets, but to under- 
stand the loss mechanisms Involved each effect 
must be separated. 

The parameters of Interest for unheated flow 
tests are the thrust coefficient Cp, the dis- 
charge coefficient Cq, and the specific Impulse 
efficiency Njsp. Figures 3 and 4 display the 
analytical and experimental values of specific 
Impulse efficiency for an orifice plate as well as 
for convergent-divergent nozzles. The specific 
Impulse efficiency decreases with decreasing 
Reynolds number for all four nozzle contours and 
the orifice. In the Reynolds number range of 500 
to 9000, Nj$p for nitrogen ranged from 80 to 90 
percent, whereas for hydrogen It ranged from 70 to 
90 percent. All of the nozzles gave the same per- 
formance within experimental error. It should be 
kept In mind that the variation of nozzle area 
ratios (120:1 to 150:1) affects nozzle efficiency. 
However, It has been shown previously that the 
variation In specific Impulse efficiency with area 
ratio amounted to only 4 percent for area ratios 
between 100:1 and 200:1 at low (500) Reynolds num- 
bers. As the area ratio Increased, the efficiency 
decreased. 1 At higher Reynolds numbers (4000), 
the variation In Nj$p for the nozzles was only 
2 percent. Although the nozzles performed the 
same within experimental error (5 percent). It 
should be pointed out that a 5-percent Increase 
In Njsp would give a 10-percent Increase In 
thrust power. Thus, a lower Input electric power 
would be required for a given target specific 
Impulse, or an Increased thrust level could be 
used. 

The difference In Nxsp for the two gases Is 
more pronounced than the difference In Ni$p for 
the various nozzle contours. Below Reynolds num- 
bers of 4000, nitrogen shows a higher than 

hydrogen, whereas at high Reynolds numbers Njsp 
Is approximately 0.88 for both gases. The maximum 
difference at a Reynolds number of 500 Is between 
5 and 12 percentage points depending on nozzle 
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configuration. At a Reynolds number of 500 f the 
greatest disparity between the specific Impulse 
efficiency of hydrogen and nitrogen was observed 
with the trumpet nozzle. 

Murch's data 1 for heated nitrogen and hydrogen 
exhibited an Njsp difference of 6 percent, which 
Is about the same observed In Figs. 3 and 4 with a 
conical nozzle. The viscous flow model of Edwards 
and Jansson 2 Indicates very little difference In 
N I$P f° r hydrogen and ammonia for Reynolds numbers 
between 200 and 10 000, Indicating a need for more 
sophisticated flow models. Comparing the different 
specific Impulse efficiencies for the orifice shows 
that, within experimental error, the performance of 
nitrogen and hydrogen Is the same. At a Reynolds 
number of 500, Nj$p for nitrogen Is 54 percent 
and for hydrogen, 57 percent. At Reynolds numbers 
approaching 8000, the specific Impulse efficiency 
for nitrogen Is 66 percent and for hydrogen, 68 
percent. Significant variation In Njsp between 
nitrogen and hydrogen only occurs when the gases 
are expanded out of the nozzle. Evidently the 
hydrogen nozzles at low Reynolds numbers have a 
thicker boundary layer than nitrogen nozzles, or 
(possibly) there Is condensation of the nitrogen 
as It expands out of the nozzle. 

Hydrogen does have a negative Joule-Thompson 
coefficient (heats upon expansion) In this 
temperature-pressure range as compared to nitrogen, 
but the effect Is too small to account for any 
N ISP variations. It Is possible that a Reynolds 
number based on the physical throat diameter Is not 
the appropriate correlation parameter. 

Figures 5 and 6 compare discharge coefficients 
Cq for both gases and all four nozzles. As can be 
seen, the bell and modified trumpet nozzles have 
higher discharge coefficients than the conical and 
trumpet nozzles (throat measurements are accurate 
to within 3 percent). The difference Is approxi- 
mately 4 percent at a Reynolds number of 8000 and 
8 percent at a Reynolds number of 500. Since the 
discharge coefficient Is a measure of the mass-flow 
rate decrease due to throat boundary-layer viscous 
effects, It appears as though the boundary layer 
Is thinner for the bell and modified trumpet noz- 
zles. For a given nozzle, the discharge coeffi- 
cients produced by hydrogen and nitrogen vary by 
no more than 2 or 3 percent. 

Back has found that the discharge coefficient 
depends essentially on the ratio of nozzle-throat 
radius of curvature to nozzle-throat radius 
r c /rth and 1s Independent of Inlet configuration 
(convergent contour). As the ratio r c /rth 
decreases, Cq decreases. 10 * 11 Although these 
results are for higher Reynolds number ranges, the 
same trends are seen at the lower Reynolds numbers, 
with the bell and modified trumpet nozzles having 
the highest ratio of r c /rth and, correspondingly, 
the highest discharge coefficients. Massler 12 has 
found no effect of plenum length (growth of bound- 
ary layer In plenum) on the discharge coefficient. 
Rothe 1 ^ has found that, for pressure ratios (gas 
pressure at nozzle Inlet to test cell pressure) 
greater than approximately 350 for a Reynolds num- 
ber of 500, the test cell pressure Influencing the 
Internal flow through the thick subsonic boundary 
layer has no effect on the discharge coefficient. 

All tests were run with this pressure ratio greater 
than TO 4 . 


All four nozzles had the same Inlet radius of 
curvature (1 mm), but their outlet radii varied 
from 0.3 to 12 mm, with the bell nozzle having the 
smallest value. The nozzles had throat length-to- 
dlameter ratios between 0.4 and 1.1, with the modi- 
fied trumpet having the smallest value and thus 
more closely approximating a sharp-edged orifice 
Inlet. The orifice, which had a discharge coeffi- 
cient equal to or greater than the best nozzle, 
had no Inlet or outlet radius of curvature 
(straight-through orifice) and had a throat length- 
to-dlameter ratio of 0.6 . It Is Interesting to 
note that the orifice was designed with no conver- 
gent section and Its discharge coefficient Is 
approximately the same as, or better than, the 
nozzles. An upstream convergent section Is not 
required to provide a high discharge coefficient 
for a choked orifice. Rae 14 predicted this fact 
analytically by using the slender channel equations 
with slip boundary conditions at the wall. It 
appears that Back's findings of the discharge 
coefficient depending on r c /r th do extend 
Into the low Reynolds number range. 10 * 11 It can 
be assumed that the throat boundary layer Is 
affected by r c /r^ f where a smaller Inlet radius 
(for a given throat radius) causes the formation 
of a larger throat boundary layer. As the throat 
boundary layer Increases, the ratio of mass flow 
through the Invlscld core to the mass flow through 
the boundary layer decreases causing a decrease In 
the average throat exit-gas velocity. 

Table 2 summarizes the basic parameters Cq, 
Nisp. and Cp for all nozzles and the orifice plate 
at a throat Reynolds number of 1000. Nozzle 
discharge coefficients are clearly Independent of 
the gases tested. At a Reynolds number of 1000, 
the highest nozzle thrust coefficient obtained was 
1.33 and the lowest was 1.14. Almost all thrust 
coefficients measured at this Reynolds number were 
less than the thrust coefficient calculated for 
Isentroplc flow through a choked orifice (Cp = 

1.27). This apparently Indicates that all nozzles 
were filled with a large boundary layer. The 
values of Cp for the orifice plate at a Reynolds 
number of 1000, 0.94 to 0.97, Indicate a signifi- 
cant departure from Invlscld flow. With nitrogen 
propellant, the variation In the nozzle thrust 
coefficient for the four different nozzle configur- 
ations was within the experimental uncertainty 
(±2.5 percent) . 

The results of the Two-Dimensional Kinetics 
(TDK) nozzle analysis can be seen In Figs. 3 and 4. 
Numerous problems were encountered during attempts 
to analyze low Reynolds number nozzles with a pro- 
gram that was written for high-thrust, high 
Reynolds number nozzles. The most serious problem 
appeared to be that the analysis would not always 
converge to a Mach number of unity at the throat. 
When running properly, the discharge coefficients 
calculated were 0.95 or greater, whereas the 
experimental data showed that the actual Cq 
varied from 0.83 to 1.0. The calculated discharge 
coefficient departed from experimental results at 
Reynolds numbers below approximately 4000. It can 
be assumed that, at this point, the TDK method of 
subtracting boundary-layer effects from the Invls- 
cld core Is nolonger valid. 

For the conical nozzle, the specific Impulse 
efficiency was predicted fairly well to a Reynolds 
number of 3000, after which substantial departure 
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from the experimental data resulted. It should be 
noted that this Is the approximate location where 
the calculated discharge coefficient began to 
deviate from the experimental results. Results 
for the bell-shaped nozzle were Inaccurate at every 
point. For the single point that was able to be 
run for the modified trumpet nozzle, the results 
were as accurate as In the conical case. No anal- 
ysis was possible for the trumpet nozzle as the 
throat Mach number would not converge to unity. 

The program did predict large displacement 
thicknesses at the nozzle exit (Fig. 7). The dis- 
placement thickness Is the distance by which the 
solid surface would have to be displaced to main- 
tain the same mass-flow rate as that of a friction- 
less flow. At a Reynolds number of 4000, the 
displacement thickness was about 40 percent of the 
exit plane area for the conical nozzle, 37 percent 
for the bell, and 67 percent for the modified trum- 
pet. The values of thrust coefficient obtained 
for the conical, bell, and modified trumpet nozzles 
at a Reynolds number of 4000 were 1.41, 1.48, and 
1.48, respectively. The calculated values of Cp 
are all large enough to Indicate the existence of 
an 1 nvl sc 1 d core, but the experimental data are 
not consistent with the TDK results. The TDK code 
needs refinements In calculating the sonic condi- 
tions at the throat for the low Reynolds numbers. 
The method of subtracting the boundary-layer vis- 
cous effects from the Invlscld core performance 
must also be refined. 

Concluding Remarks 

Four nozzles (a conical, bell, trumpet, and 
modified trumpet) and a sharp-edged orifice were 
evaluated over a Reynolds number range of 500 to 
9000 with unheated nitrogen and hydrogen. The 
boundary-layer viscous losses were examined In 
both the divergent section and throat section of 
the nozzles. The nozzles showed significant 
decreases In specific Impulse efficiency with 
decreasing Reynolds number. We conclude that 
changes In the divergent contour do not effect the 
viscous and divergent losses to any appreciable 
extent. At Reynolds numbers less than 1000, all 
four nozzles were probably filled with a large 
boundary layer. 

There appears to be an effect of gas species 
on the viscous and divergent losses with Reynolds 
number. At a Reynolds number of 500, nitrogen 
performs approximately 5 to 12 percent better than 
hydrogen, depending on nozzle configuration. At 
the higher Reynolds numbers (>4000), the difference 
Is within 3 percent. The difference In specific 
Impulse efficiency does not appear when testing the 
orifice. This difference Is related to the expan- 
sion process taking place In the divergent portion 
of the nozzles. 

The discharge coefficient decreases with 
Reynolds number In the same manner as the specific 
Impulse efficiency, but the gas species has very 
little effect on the discharge coefficient. The 
bell and modified trumpet nozzles, with higher 
ratios of Inlet radius of curvature to throat 
radius, have 4 to 8 percent higher discharge coef- 
ficients than the cone or trumpet nozzles. 

A comparison of the experimental results to 
the Two-Dimensional Kinetics (TDK) nozzle analysis 
predicted results Indicates that better predictions 


for the discharge coefficient are needed along 
with refinements In calculating the sonic condi- 
tions In the nozzle throat for Reynolds numbers 
below approximately 4000. The Incorporation of 
viscous effects by subtracting the boundary-layer 
viscous losses from the Invlscld core performance 
calculations Is Insufficient to account for the 
large viscous losses at the low Reynolds numbers. 

A more adequate method of accounting for viscous 
losses Is needed. 

Although the divergent section geometry has 
little effect on the specific Impulse efficiency, 
a 5-percent Increase would provide a 10-percent 
Increase In thrust power per thruster. This would 
amount to a substanlal payback for a satellite 
system composed of numerous thrusters. 

The knee of the specific Impulse efficiency 
curve falls In the region that Is representative 
of the Reynolds number regime In which high- 
performance electrothermal thrusters operate. 

This region Is Important for two reasons. Satel- 
lite systems operate In a blowdown mode where the 
propellant flow rate decreases with use. Conse- 
quently, as the satellite ages Its performance will 
decrease (following the specific Impulse effi- 
ciency curve). There Is also a trade off (for a 
given flow rate) between specific Impulse effi- 
ciency and thruster chamber high-temperature 
creep. Large nozzles (large throat diameters) 
decrease the thruster operating pressure, thus 
decreasing the creep problems of high-temperature 
materials. However, there Is a corresponding 
decrease In Reynolds number and, therefore, a 
decrease In specific Impulse efficiency. 

More analytical work needs to be done In this 
Reynolds number regime to aid In evaluating the 
above tradeoffs from a systems standpoint. 
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TABLE 1. - NOZZLE PARAMETERS 


Nozzle 

Shape 

Exit half-angle 

Throat diameter, 
mm 

Area ratio 



Throat 

Exit plane 


1 

Conical 

20° 

20° 

0.653 

120:1 

2 

Bell 

35° 

o 

O 

OJ 

.711 

150:1 

3 

Trumpet 

0° 

36° 

.671 

125:1 

4 

Modified 

trumpet 

15° 

o 

GO 

OJ 

.640 

135:1 


TABLE 2. - NOZZLE CHARACTERISTICS AT A REYNOLDS NUMBER 
OF 1000 FOR NITROGEN AN0 HYDROGEN 


Nozzle 

Shape 

Discharge 
coef f Iclent, 
CD 

Specific Impulse 
efficiency, 

N ISP 

Thrust 
coef flclent , 
Cf 



Nitrogen 

Hydrogen 

Nitrogen 

Hydrogen 

Nitrogen 

Hydrogen 

1 

Conical 

0.87 

0.88 

0.80 

0.78 

1.18 

1.17 

2 

Bell 

.91 

.92 

.85 

.75 

1 .31 

1.17 

3 

Trumpet 

.88 

.86 

.83 

.78 

1.24 

1.14 

4 

Modified 

trumpet 

.93 

.93 

.84 

.76 

1.33 

1.20 

- 

Orifice 

plate 

.98 

.94 

.58 

.59 

.97 j 

.94 
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FIG. 1. - NOZZLE GEOMETRIES. ALL DIMENSIONS ARE IN MILLIMETERS. 



CALIBRATION 

WEIGHTS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


r 


QUICK-CONNECT 



/' FLANGE ASSEMBLY 

/ 


WINDSHIELD 


L LI NEAR 
VARIABLE 
DISPLACEMENT 
TRANSDUCER 
(LVDT) 


TEMPERATURE ■ 


ELECTROMAGNETIC , 
DAMPER AND / 
MOVABLE MAGNETS J 


I 


NOZZLE 


O-RING SEAL 


-THRUSTER 
MOUNTING PLATE 


V FLEXURE PLATES 


1 


(A) SCHEMATIC OF NOZZLE ASSEMBLY AND THRUST STAND. 



(B) THRUST STAND. 

FIG. 2. - NOZZLE ASSEMBLY AND THRUST STAND. 
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THROAT REYNOLDS NUMBER 


FIG. 3. - SPECIFIC IMPULSE EFFICIENCY AS FUNCTION OF 
THROAT REYNOLDS NUMBER FOR HYDROGEN. EXPERIMENTAL 
DATA AND TWO-DIMENSIONAL KINETICS (TDK) NOZZLE 
ANALYSIS. 



THROAT REYNOLDS NUMBER 


FIG. 4. - SPECIFIC IMPULSE EFFICIENCY AS FUNCTION OF 
THROAT REYNOLDS NUMBER FOR NITROGEN. EXPERIMENTAL 
DATA AND TWO-DIMENSIONAL KINETICS (TDK) NOZZLE 
ANALYSIS. 
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THROAT REYNOLDS NUMBER 


FIG. 5. - DISCHARGE COEFFICIENT AS FUNCTION OF THROAT 
REYNOLDS NUMBER FOR HYDROGEN. 



THROAT REYNOLDS NUMBER 


FIG. 6. - DISCHARGE COEFFICIENT AS FUNCTION OF THROAT 
REYNOLDS NUMBER FOR NITROGEN. 


TO 



2.17 



THROAT REYNOLDS NUMBER 

FIG. 7. - TWO-DIMENSIONAL KINETICS (TDK) NOZZLE ANALYSIS 
PREDICTED DISPLACEMENT THICKNESS AS FUNCTION OF REYNOLDS 
NUMBER AT EXIT PLANE FOR NITROGEN AND HYDROGEN. 
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